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SUMMARY 


A rigid flap-lag blade analysis has been developed to simulate a rotor in a wind tunnel under- 
going an emergency power shutdown. Results show that for a rotor at a nonzero shaft tilt angle 
undergoing an emergency power shutdown, the oscillatory lag response is divergent. The mean lag 
response is large when tested at high collective pitch angles. Reducing the collective pitch during the 
emergency shutdown reduces the steady lag response. Increasing the rotor shaft tilt angle increases 
the oscillatory lag response component. The blade lag response obtained by incorporating a nonlin- 
ear lag damper model indicates that in this case the equivalent linear viscous damping is lower than 
originally expected. Simulation results indicate that large oscillatory lag motions can be suppressed 
if the rotor shaft is returned to the fully vertical position during the emergency power shutdown. 


INTRODUCTION 


This report describes the rigid-blade analysis of a hovering rotor while it is undergoing an 
emergency shutdown in rotor power. The results of this report support the safety analysis of the 
Sikorsky S-76 rotor tested in Ames Research Center’s 80- by 120-Foot Wind Tunnel. 

During a wind tunnel test of the Sikorsky S-76 rotor, the rotor shaft can tilt as high as 30°. This 
test condition is necessary to minimize ground effects, and it poses several safety concerns. If an 
emergency shutdown in rotor power occurs at these critical testing conditions, then the blade lag 
motions can become large. Large blade lag motions may lead to two drastic consequences: 

(1) fatigue failure of structural components, and (2) ground resonance. This report analyzes the 
rotor blade motions subsequent to an emergency shutdown in rotor power. 


Ground Resonance - Normal Operating Conditions 

Most helicopter rotors tested in a wind tunnel can experience ground resonance. Ground reso- 
nance is a violent dynamic instability caused by the coupling between the inplane hub and blade 
motions, and it occurs only in a certain critical range of the rotor rotational speed. The installation of 
blade lag dampers is crucial to prevent ground resonance. For this, the damping level is typically 
sufficient only for the normal operating conditions. 


Ground Resonance - Critical Operating Conditions 

After an emergency shutdown in rotor power, the rotor speed reduces rapidly. Gravity then 
becomes a dominant load in the blade aeroelastic responses and can generate large blade lag 
motions. This is a serious problem, especially for rotors tilted at large shaft angles. 



For the Sikorsky S-76 rotor, the lag damper performance is significantly reduced when the blade 
lag motions exceed a specific value, and if this happens at the critical rotor speed range, then the 
rotor becomes susceptible to ground resonance. 
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Time constant for rotor speed decay, sec 
Lift curve slope 

Time constant for dynamic inflow model 
Blade chord, ft 

Equivalent viscous damping coefficient, lb-ft/sec 

Drag coefficient, constant term 

Drag coefficient, quadratic term 

Hinge offsets for both flap and lag, ft 

Gravitational constant 

Flapping moment of inertia, slug-ft 2 

Flap hinge spring constant, lb-ft 

Lag hinge spring constant, lb-ft 

Pitch-flap coupling (pitch up-flap up is positive) 

Pitch-lag coupling (pitch up-lag back is positive) 

Lag damper moment, lb-ft 

Control rate, deg/sec 

Blade radius, ft 

Blade first moment of inertia, slug-ft 
Time, sec 
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vj Rotor inflow, ft/sec 

Xh Longitudinal hub displacement (positive aft) 

<Xs Rotor shaft tilt angle (positive aft), deg 

P Blade flap (positive up), rad 

Pp Blade precone 

C, Blade lag (positive aft), rad 

Blade pre-lag 

Q Rotor rotational speed 

0.75 Blade collective pitch at 0.75 blade radius 

0lc Lateral cyclic pitch 

01s Longitudinal cyclic pitch 

0 tw Blade linear pre-twist 

a Rotor solidity 

( ) Time derivative 


ANALYTICAL FORMULATION 


The analysis is developed for a rotor experiencing a temporal variation in its rotational speed. 
The rotor blade is modeled as a rigid beam undergoing flap (out-of-plane) and lag (inplane) degrees 
of freedom. The modeling is valid for rotor blades with linear pre-twist, pre-cone, pre-lag, and kine- 
matic pitch-flap and pitch-lag couplings. Linear quasi-steady aerodynamics is used to represent the 
blade airloads. Dynamic inflow model based on momentum theory is used to calculate the rotor- 
induced inflow. The effects of gravity and longitudinal shaft motions are included. Hub flexibility 
and forward speed effects are not included. 

To simulate the rotor power shutdown, the rotor rotational speed is assumed to undergo an expo- 
nential decay. The time constant value for the decay rate is taken from the results of the Tunnel 
Utilization Trainer with Operating Rotor (TUTOR) analysis (ref. 1). This simplification captures the 
effects of rotor shaft torque dynamics along with motor engine and gearbox dynamics of the drive 
system for the Ames Rotor Test Apparatus. With this assumption, the rotational speed becomes 
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( 1 ) 


Q(t) = Q 0 e -t/A i 

The governing equations for a rigid flap-lag blade are derived by integrating individual compo- 
nents of the blade section moment about the flap and lag hinges. For the Sikorsky S-76 rotor blade, 
the flap and lag hinges coincide, and therefore the offset value, e, applies for both hinges. These 
equations are then normalized with respect to lb, the blade flap moment of inertia. Note that these 
equations are not normalized with Ib£2 2 as is the standard practice due to the temporal variation of Q. 
The resulting equation for flap is 
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and that for lag is 
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where 


g x = - g sina s /R 
g y = g cosa s /R 
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xh = 0.4745 a s sin'F 


01 = 0.75 + 01c cos'? + Oissin'F - O.750 tw + KpP + 

= - 0.01 25^ - 0.00375P + 0.0002500.75 + 0.07875 

vp = 1 + — f m(r-e) 2 dr 
lb 

Vr = — f m(r - e) 2 dr 
•'e 

f R 2 

I b = I m(r-e) dr 
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S b = m(r-e)dr 

Je 

°Jpo = k p 1 h 

o = k c 7I b 

S b = R s b /l b 

- do / a 
d2 = d2 / a 
A. = Vj/R 

'p = f a dt 
Jo 

Note that the pitch-lag coupling term is taken from reference 2 and depends on the blade flap 
and lag motions as well as the collective pitch angle. The effects of blade flap and lag stops are also 
included. Due to the lack of information about the stop dynamics, it is assumed that the blade flap or 
lag motion stops after the blade hits the corresponding stop. 


The Sikorsky S-76 rotor employs nonlinear lag dampers. In the present analysis, the effects of 
both a linear viscous damper model and a nonlinear damper model are studied. For the linear viscous 
damper. 


£ 


For the nonlinear damper (ref. 3) 
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M£ = 


19482 £ 


rl 


378 + 454 




for 


for 


|5.84 £ < 1 


5.84 0 > 1 


The quantity 5.84 £ denotes the lag damper velocity in in./sec. From the preceding expression for 
the nonlinear lag damper model, note that the lag damper performance is drastically reduced when 
the damper velocity exceeds 1 in./sec. Also let 

= Mj- / Ijj 

c^q./ib 

The rotor-induced inflow is obtained from momentum theory, and the dynamic effects are 
included with a first-order time lag. The dynamic inflow equation is 


1*1 ]C y 
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The value of a i is assumed to be 2.0. 

Equations (2)-(4) are combined and expressed in a first-order form, and the resulting first-order 
differential equation is numerically integrated in time using the Gear method (ref. 4). 


RESULTS AND DISCUSSIONS 


Results are obtained using the properties of the Sikorsky S-76 rotor (ref. 2). The blade and rotor 
properties used in the simulation are presented in table 1 . Figure 1 shows the blade flap and lag fre- 
quencies as functions of rotor rotational speed. Due to the matched-stiffness design, the nonrotating 
flap and lag mode frequencies are equal. 

To initiate the time integration, initial conditions have to be prescribed for the rotor degrees of 
freedom and inflow. For this, a preliminary simulation is performed at a constant rotor speed, and the 
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initial conditions for the rotor states (blade flap and lag displacements and velocities, and the rotor 
inflow) are prescribed arbitrarily. The steady state results yield the initial conditions for the actual 
simulation when the rotor speed decays. 

Figure 2 shows the prescribed decay value of the rotor speed as the power is cut off at zero sec, 
and this rotor speed variation is used in all subsequent simulations. 

The first two sets of results represent the shaft-fixed conditions. They are based on the lag 
damper modeling: (1) linear viscous model, and (2) nonlinear model. The third set of results shows 
the effects of shaft motions. 


Linear Viscous Lag Damper Model 

Effect of rotor shaft tilt- Figures 3 and 4 show the effects of shaft tilt on blade responses under 
a power shutdown condition. The results are presented at three shaft tilt angles: 0, -15, and -30°. 

The other test parameters and lag damper values are presented in table 2. Figure 3 shows the blade 
flap responses. Clearly, the blade flaps move downward as rotor speed decreases. After about 30 sec, 
the blade hits the flap stop at -6° for the three cases shown. Except for the higher harmonics which 
increase with the higher shaft tilt angle, all the flap responses behave similarly. The higher harmonic 
flapping motions are due to the Coriolis coupling between the flap and lag motions. 

The lag responses are shown in figure 4. The results are presented about the pre-lag angle of 7°. 
The total blade lag is that shown, as in figure 4 and subsequent figures as well, with an additional 7° 
of pre-lag. The mean lag (the steady components) is slightly reduced at the higher shaft tilt angles. 
The oscillatory lag (the higher harmonic components) is moderate at a shaft tilt angle of 15° and 
increases at a shaft tilt angle of 30°. Note that the oscillatory lag motion is absent at zero shaft tilt 
angle. 

The blade lag response can be considered as a forced response since its general behaviors are 
dependent on the external forcing. Furthermore, the above results imply that the total lag response 
can be separated into two components: mean and oscillatory components. These two facts suggest 
that (1) the rotor speed variations govern primarily the mean lag response behavior, and (2) both the 
rotor speed variation and the component of gravity acting in the rotor inplane direction govern the 
oscillatory lag response. 

Effect of collective pitch- The effect of collective pitch on blade responses under a power-loss 
condition are shown next. For these results, four different values of the blade collective pitch are 
used, ranging from 5 to 18°. The collective pitch remains constant during the power shutdown pro- 
cess. Other parameters used in this simulation are presented in table 2. The blade flap responses 
shown in figure 5 indicate that different values of collective pitch have different initial blade 
responses. However, the decay behaviors for these cases are quite similar, and for all four cases 
shown, the blade hits the flap stop at roughly 30 sec after the power shutdown starts. 

Figure 6 shows the lag responses. Even though the steady components of lag response are quite 
different for different values of collective pitch, they all have a similar oscillatory behavior. When 
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the blade collective pitch is above 15°, the blade hits the lag stop at 10 ° (a value measured about a 
pre-lag angle of 7°) about 32 sec after the power shutdown. From the safety point of view, this is not 
a major concern since this happens when the rotor speed is low, at about 30 rpm (fig. 2). 

Effect of collective pitch reduction rate- In the following simulations, the collective pitch 
varies during the power shutdown. The control schedule is 


0 75 (t) = 075 o -r*t ( 6 ) 

and (0 75 ) . is 2 °. In the above expression, r is the control rate, t denotes time in sec, and ©75 q is 
the initial coYlective pitch value. The following set of results is shown in figures 7 and 8 for four dif- 
ferent values of r, ranging from 0 to 20 deg/sec. The value of r for the Ames Rotor Test Apparatus is 
2 deg/sec. Table 2 presents the values of other parameters used. 

The blade flap responses shown in figure 7 indicate that higher values of r result in a shorter 
decay time. However, the flap responses are essentially the same for the values of r higher than 
2.0 deg/sec, except for the short transient behaviors. The similarity in steady state behaviors are 
expected since the pitch schedule with the higher pitch rate (r > 2 deg/sec) are equal after 5 sec. 

For the lag response results shown in figure 8 , reducing the collective p : tch after a power loss is 
beneficial because this leads to a smaller steady lag motion and can prevent tne blade from hitting 
the lag stop. The oscillatory components are not affected by the control reduction rate. As in the case 
with flap response, there is a limiting value to the control rate effectiveness. 

Effect of linear viscous damping- The blade responses with different damping values are 
shown in figures 9 and 10. Other parameters used are also presented in table 2. As shown in figure 9, 
the flap responses are quite insensitive to the blade lag damper values. 

The lag responses shown in figure 10, on the other hand, are quite sensitive to the lag damper 
values. The steady lag responses decrease slightly, up to 2° maximum, with a larger damping value. 
The oscillatory lag amplitudes also decrease with an increase in lag damping value. 


Nonlinear Lag Damper Model 

The study on the effects of rotor shaft tilt and collective pitch on the blade responses under a 
power shutdown are repeated with the nonlinear lag damper model. The primary test values are 
shown in table 2. For these results, the flap responses are not shown since, as expected they are 
found to be similar to those obtained with the linear lag damper model. 

Effect of rotor shaft tilt- Figure 1 1 shows the effects of rotor shaft tilt on the blade lag 
responses. Even though the general behavior is similar to the results with linear lag damper model 
(fig. 4), the actual responses are quite different. At zero shaft tilt angle, the lag response still does not 
exhibit any oscillatory behavior. However, the blade hits the lag stop at 35 sec compared to 42 sec 
for the same result with the linear lag damper model. At the nonzero shaft tilt angles, the components 
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of oscillatory lag response are divergent. The oscillatory magnitudes are almost twice those obtained 
with the linear lag damper model, while the oscillating frequencies remain the same. 

Effects of collective pitch- The effect of collective pitch on the blade lag responses with the 
nonlinear lag damper model are presented in figure 12. Comparing these results to those of figure 6, 
it is seen that there are two consequences of the nonlinear lag damper model: (1) the blade hits the 
lag stop earlier, and (2) the oscillatory lag responses are almost twice that of the results shown in 
figure 6. 

Linear vs. Nonlinear Lag Damper Models- The blade lag responses with a linear lag damper 
model are compared with that with the nonlinear lag damper model in figure 13. For the linear lag 
damper results, two values of the damping coefficient are used: 2000 and 6000 ft-lb-sec. From fig- 
ure 13, it can be seen that the lag response using the linear lag damper model with a damping coeffi- 
cient of 2000 ft-lb-sec is identical to that using the nonlinear lag damper model. This suggests that an 
equivalent linear viscous damper model with an appropriate damping coefficient can represent a 
nonlinear damper. The results of figure 13 also indicate that the linear lag damper model with a 
damping coefficient of 6000 ft-lb-sec is an overestimate of the equivalent viscous damping. 


Effects of Rotor Shaft Motions 

All the previous results are obtained with the rotor shaft fixed during the emergency power 
shutdown. For the results of this section, the shaft tilt angle is varied by the expression 

a s (t) = a Sj0 +a s *t (7) 

where a Si o is the initial shaft tilt angle and a s is the shaft tilt rate. The shaft tilt rate is assumed to be 
zero once the shaft is brought to the fully upright position, i.e., a s = 0°. The shaft tilt rate for the 
Ames Rotor Test Apparatus is assumed to be 1 deg/sec. This simulation uses the nonlinear lag 
damper model, and table 2 presents the values of other parameters used. The flap responses are 
identical to those of figure 3 and are not shown. 

Figure 14 shows the lag responses obtained at different initial shaft tilt angles. The results in this 
figure show that the oscillatory lag responses are not divergent and are much smaller than those 
shown in figure 1 1 , corresponding to the nonzero shaft fixed conditions. In fact, the oscillatory lag 
components all vanish before the blade hits the lag stop. From the safety point of view, these results 
indicate that large oscillatory lag motions can be avoided if the rotor shaft is returned to the fully 
vertical position during the emergency power shutdown. 
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CONCLUSIONS 


An analysis based on a rigid flap-lag blade model has been developed and used to predict the 
aeroelastic responses of a hovering rotor undergoing an emergency shutdown in rotor power. A 
simulation study is carried out for different rotor properties and operating conditions. 

The results indicate that for a rotor tested at a nonzero shaft tilt angle, the oscillatory components 
of lag response are divergent and can be reduced with an increase in the lag damping. Reducing the 
collective pitch during the power shutdown reduces the steady component of the lag response. 
Results with a nonlinear lag damper model suggest that an equivalent linear damper with a damping 
coefficient of 2000 ft-lb-sec is appropriate to represent the Sikorsky S-76 rotor damper. If the rotor 
shaft is returned to the fully vertical position during the emergency power shutdown, then the large 
oscillatory lag motions can be avoided. 
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Table 1. Blade and rotor properties used 
in the simulation model. 


A = 13.16 sec 
a = 5.73 
aj =2.0 sec 
e = 0.833 ft 
d 0 = 0.0087 
d 2 = 0.1719 
c= 1.29 ft 
I b = 408 slug-ft 2 
kp = 1 192 ft-lb 
k^= 1192 ft-lb 
Kp = -0.305 
R = 22.0 ft 
Sb = 30 slug-ft 


Pmax — 21 deg 
Pmin = “6 deg 
Pp = 0.0 deg 
y= 10.09 
np = 1.0296 
V£ = 0.2470 
a = 0.0748 
e tw = -10 deg 
01 c = 0.00 deg 
0] s = 0.00 deg 
Cmax =17 deg 
Cmin = “5 deg 
= 7 deg 


Table 2. Values of test parameters used in the simulations. 


Figures 

a, deg 

075 . deg 

r, deg/sec 

C^, ft-lb-sec 

3,4,11 

Varies 

10 

0.0 

6000* 

5, 6, 12 

-15 

Varies 

0.0 

6000* 

7,8 

-15 

10 

Varies 

6000* 

9, 10 

-15 

10 

0.0 

Varies 

13 

-15 

10 

0.0 


14 

Varies 

10 

0.0 



*For linear lag damper model only. 
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Figure 3. Effects of shaft tilt on flap response of rotor undergoing an emergency shutdown in rotor 
power. 
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Figure 4. Effects of shaft tilt on lag response of rotor undergoing an emergency shutdown in rotor 
power. 
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Figure 5. Effects of collective pitch on flap response of rotor undergoing an emergency shutdown in 
rotor power. 
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Figure 6. Effects of collective pitch on lag response of rotor undergoing an emergency shutdown in 
rotor power. 
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Figure 7. Effects of collective pitch rate on flap response of rotor undergoing an emergency 
shutdown in rotor power. 
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Figure 8. Effects of collective pitch rate on lag response of rotor undergoing an emergency 
shutdown in rotor power. 
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Figure 9. Effects of lag damping amplitude on flap response of rotor undergoing an emergency 
shutdown in rotor power. 
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Blade lag, deg 



Figure 10. Effects of lag damping amplitude on lag response of rotor undergoing an emergency 
shutdown in rotor power. 
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Figure 1 1 . Effects of shaft tilt on lag response of rotor undergoing an emergency shutdown in rotor 
power-nonlinear lag damper model. 
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Figure 12. Effects of collective pitch on lag response of rotor undergoing an emergency shutdown in 
rotor power-nonlinear lag damper model. 
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Figure 13. Effects of lag damper model on lag response of rotor undergoing an emergency shutdown 
in rotor power. 
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Figure 14. Effects of shaft motion on lag response of rotor undergoing an emergency shutdown in 
rotor power. 
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